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Extreme Electron-Phonon Coupling in Boron-based Layered Superconductors 
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The phonon-mode decomposition of the electron-phonon coupling in the MgB2-like system 
Lii-zBC is explored using first principles calculations. It is found that the high temperature su- 
perconductivity of such systems results from extremely strong coupling to only ~2% of the phonon 
modes. Novel characteristics of Eig branches include (1) "mode A" values of 25 and greater com- 
pared to a mean of ~ 0.4 for other modes, (2) a precipitous Kohn anomaly, and (3) Eig phonon 
linewidths within a factor of ~2 of the frequency itself, indicating impending breakdown of linear 
electron-phonon theory. This behavior in borne out by recent inelastic x-ray scattering studies of 
MgB 2 by Shukla et al. 



Superconductivity near T c w 40 K in MgB2 [1] 
has necessitated a re-evaluation of our understanding 
of phonon-coupled superconductivity, and illustrated 
vividly that there has been a substantial void in our 
conceptual understanding of this "conventional" mech- 
anism of coupling, long thought to require d electrons, 
high density of states, and high symmetry. The gen- 
eral theory of electron-phonon (EP) coupled supercon- 
ductivity has long been available, [2,3] and some essen- 
tial understanding of the mechanism in MgB2 is now in 
place. [4-10] The particular structure and chemistry of 
MgB2 creates holes in the B 2pa bands, and these holes 
are coupled very strongly to the B-B bond stretching 
vibrational modes. The predictions of band theory have 
been verified by angle-resolved photoemission [11] and de 
Haas - van Alphen measurements. [12-14] One new facet 
is an unusually strong anisotropy in the EP coupling, 
such that the Fermi surface (FS) averaged EP coupling 
strength A and Eliashberg function a 2 F(w) do not give 
a precise account of T c or of thermal and spectroscopic 
data. It seems however that much can be accounted for 
by generalizing to a two-band model [9,10] of strongly 
coupled a holes and weakly coupled n electrons, with 
A = A CT + A,,- ~ 0.6-0.8. 

In spite of success in accounting for many observations, 
this linear EP theory (with some ad hoc correction for 
anharmonicity) seems at best incomplete. The enormous 
linewidth 7 of the Ei g mode, [15-18] ascribed to strong 
EP coupling plus anharmonicity, has not been accounted 
for quantitatively, and the predicted change in frequency 
and 7 at T c [9] is not observed. Non-adiabatic processes 
arising from the small filling of the a bands, and low 
Fermi energy E^, have been proposed. [9,19,20] Boeri et 
al. provided needed clarification by demonstrating that 
E2 g anharmonicity arises in MgB2 from non-adiabatic 
effects, [21] due to the huge deformation potential. [4] 

A second compound in the "MgB2 class" that has 
been proposed recently is useful in clarifying these issues. 
LiBC is isostructural [22] and isovalent with MgB2, hav- 
ing graphene-like B-C layers that are even more strongly 



bonded than are the B layers of MgB2 . The stoichiomet- 
ric compound is semiconducting (due to the B-C differ- 
ence) , but Li deficiency that makes it a hole-doped metal 
has been reported. [22] Such hole doping introduces car- 
riers into the B-C a bands, as occurs in MgB2, and near 
rigid-band behavior has been demonstrated, making it a 
promising system for probing the carrier density depen- 
dence of EP coupling in MgB2-class superconductors 

In this paper we focus on the decomposition of A 
into contributions ("mode A's") from each phonon Qv 
rather than its FS decomposition, and find startlingly 
high values \q v w 25 for Ei g modes. First principles 
calculations of a 2 F for Lio.75BC gives A=0.74, making 
Lii_a;BC at this doping level quantitatively similar to 
MgB2. The insulating end member of the Lii-^BC sys- 
tem makes it pedagogically an ideal material to illustrate 
the novelty that MgB2 has introduced into the physics of 
electron-phonon coupled superconductivity. Lii-^BC is 
"MgB2 with stronger coupling and variable carrier con- 
centration," and because its a FSs are (as in MgB2) 
very close to cylindrical, the underlying physical pro- 
cesses can be modelled simply yet realistically. The es- 
sential revelations are: (1) the high temperature super- 
conductivity arises from exceedingly strong coupling of 
only a small fraction (~2%, depending on doping) of the 
phonon modes, (2) a colossal sharp Kohn anomaly occurs 
- extreme renormalization by EP coupling upon doping 
from hard bond-stretching modes, and (3) the phonon 
linewidth due to decay into a band electron-hole pairs 
becomes comparable to the linewidth itself, making the 
bond-stretching modes poorly defined (at best), yet not 
unstable (due to their intrinsic hard nature). Just this 
strength of coupling and relative linewidths have been 
reported for MgB2 by Shukla et al., obtained from in- 
elastic x-ray scattering spectrosciopy. [23] Moreover, (4) 
these effects occur continuously with doping and will al- 
low a unique window into the non-adiabatic regime where 
the Fermi energy is comparable to the phonon frequency. 
These features apply directly to Mgi_ a; Al 2: B2 as well. 

A description of the band structure of Lii-^BC was 



presented earlier. [24] The FSs for x=0.25 pictured in 
Fig. 1 , consist of four B 2p - C 2p a band cylinders with 
very small k z dispersion, and a n band FS. The a cylin- 
ders are similar to those in MgB2, with four (instead of 
the two in MgB2) arising from the doubling of the unit 
cell along c due to alternate stacking of B-C layers. The 
phonon energies and EP matrix elements have been ob- 
tained from linear response theory, [25] as implemented 
in Savrasov's full-potential linear muffin-tin orbital code. 
[26,27] Because of the emphasis here on specific Fermi 
surface effects, a dense grid of Q points was chosen (a 
16,16,4 grid giving 90 Q points in the irreducible wedge). 
For k-space integration a finer 32,32,8 grid was used, to- 
gether with an adaptive tetrahedron integration scheme. 
[26] The code was used previously for MgB2 by Ander- 
sen's group, [6] where the need for careful zone sampling 
was also emphasized. 

For the semiconducting phase, the E2 9 phonons (Fig. 2) 
are very hard, lying at 150 meV at the T point and not 
dispersing greatly. This high bond-stretching frequency, 
twice as high as in MgB2, is due to the much stronger 
B-C bond, 9% shorter and having a 40% larger defor- 
mation potential for the a bands. [24] The spectrum is 
indicative of a stable, strongly bonded material. For x = 
0.25, where it has become a metal (Fig. 2), there is rather 
little difference in the spectra, except for Qn < 2kp w 
7r/3a (Q|| is the component of the phonon wavevector Q 
lying in the plane). Within this region, however, the four 
£^2 S -related modes display a precipitous Kohn anomaly 
at 2fcp, and dip nearly dis continuously by 40%- (The 
decrease in u> 2 , which directly reflects the renormaliza- 
tion, is nearly 65%). Closely related renormalization of 
phonon modes in the Mgi-^ALjT^ system has been noted 
by Renker et al. [28] 

In Fig. 3 the phonon density of states (DOS) for x = 
and x = 0.25 are contrasted, and also the calculated EP 
spectral function a 2 F(uS) for x = 0.25 is compared with 
the DOS. For the DOS, doping moves the modes in the 
frequency range 135 - 155 meV to lower frequency; only 
from the dispersion curves of Fig. 2 is it clear where they 
end up (in the 85 meV region) because these modes are 
near the T point where phase space is small. The curve 
for a 2 F(uj) demonstrates what results: a tiny fraction of 
modes (quantified below) in the 85 meV region are ex- 
tremely strongly coupled to the a band holes, and leads 
to an inordinately strong contribution to A, whose total 
value calculated from a 2 F is 0.74. Although about half 
of the coupling arises from all other modes, it is the ex- 
tremely strong coupling of small-Q F, 2g modes to a holes 
that drives T c , just as in MgB 2 . The calculated value of 
T c using the Allen-Dynes equation [29] is 34 K assum- 
ing /j,* = 0.09 as has been used previously [10] to obtain 
numerical agreement with the observed T c for MgB2. So- 
lution of the anisotropic Eliashberg equations [10] would 
give T c a few degrees higher. 



The circular a Fermi surfaces of MgB 2 and related ma- 
terials allow an analytic treatment of EP coupling. The 
phonon energies (h=l) are given by the poles of the lat- 
tice Green's function 



O 2 
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(1) 



where the reference frequency VIq v includes all self- 
energy effects except those arising from electron scatter- 
ing within the a bands; practically speaking, they are 
the energies shown in Fig. 2 for undoped LiBC. Tak- 
ing for simplicity a single cylindrical Fermi surface with 
N(e)=m*/27r per spin per unit area, and considering the 
E 2g derived modes only, the phonon self-energy due to 
the a carriers is (t]q = Qu/21cf) 

TL'(Q,u) = -2 V|M fe , Q | 2 iWfc+Q . x 

£k+Q -Sk-uj-id 

5Re W(Q, UQ ) « -2\M 2g \ 2 N{e F )xT{Q), (2) 

X 2 L D (Q) = 6(1 - VQ ) V Q 2 ]d(riQ - !)• (3) 

Here fk = /(efe) is the Fermi occupation factor, and the 
usual adiabatic approximation has been made. A mean 
square matrix element has been extracted from the sum, 
leaving the unitless Lindhard function Xl°- [30] The neg- 
ative of Xl D {Q) nas much the form that is evident in the 
softened T^ 2g branches with Q < 2fcp in the bottom panel 
of Fig. 2. This behavior is illustrated by plotting ujq 
given by Eqs. 1-3 in the bottom panel of Fig. 4. 

The "mode A" Aq„ [31] for the E 2g modes, defined 
such that it is an intensive quantity whose average over 
the zone and over the N„ branches gives \ 2g (A = 

■Mother 

), is 

HQ) = £ s N^^MBll (4) 

where £(r]) = ?7 _1 (1 — r] 2 ) _1 / 2 6'(ry)0(l — n). Here Abz is 
the basal area of the BZ and Aps = n(2kp) 2 is a doping- 
dependent area of renormlalized E 2g phonons. Q — > 
requires somewhat more care [32] but involve negligible 
phase space for the present purposes. The important fea- 
ture, aside from the magnitude, is that these mode A's 
scale inversely with k 2 F as illustrated in Fig. 4, diverging 
for vanishing a hole doping and indicating the breakdown 
of conventional EP theory (see below) . 

For Q < 2k f, wq and are nearly constant, so Eqs. 
(1) and (2) can be solved for |M 2s | 2 and substituted into 
Eq. (4) to provide a simple expression for the mode Aq 
(Q < 2kp) in terms of (fi 2 ^ — o->2 g )/o-i 2 g - Similar ex- 
pression have been used for Nb-Mo alloys [33] and for 
Pb alloys [34] to relate A to phonon softening. There 
is a simpler way to obtain the average mode A for E 2s 



phonons. The total contribution is that obtained from 
the 85-100 meV region (see Fig. 3) of a 2 F, \ 2g =0.4. 
This arises from an average over all phonons, but only 
those in 3/40 of the BZ (Q < 2&f) and the upper four of 
18 of the branches contribute. Hence the mean is < Xq > 
= 0.4x(40/3) x (18/4) w 25. The full Q dependence of 
the mode Xq for the E 2g branches is given by Eq. (4) 
and shown in the top panel of Fig. 4. The mean mode 
A from all other phonons (equal to their contribution to 
A since it includes 98% of all modes) also is Xother ~ 0.4, 
accounting for the total A = 0.74 obtained from a 2 F. 

The phonon relative halfwidth [31] becomes jq/ivq <~ 
0.5Aq/ < X 2g >^ 0.5, an alarming result (this ratio is 
typically 10~ 2 - 10~ 3 ) that reflects the fact that this 
phonon branch is so ill-defined from extremely strong EP 
coupling that Migdal theory (and Eliashberg theory) is 
no longer justified. Even so, this rough estimate is consis- 
tent with the large observed halfwidth [15-18] of 125-175 
cm -1 for the 600 cm -1 bond-stretching mode, indicat- 
ing that the large linewidth and its strong increase with 
temperature is due to EP coupling rather than to anhar- 
monicity, and the same conclusion has been reached by 
Shukla et al. [23] 

Now we summarize and discuss some implications of 
these results. Use of Lii-^BC has allowed us to identify 
and quantify the drastic phonon softening arising from 
the ultrastrong EP coupling to E 2g modes; comparable 
phenomena occur also in MgB 2 . Mode A2 9 's w 25 and 
linewidths comparable to the frequency point to inade- 
quacies of linear EP theory for these systems. Shukla 
et al. have reported [23] measurements of specific fre- 
quencies and linewidths of MgB2 from inelastic x-ray 
scattering. For five small-Qn F> 2g phonons with frequen- 
cies in the 500-550 cm -1 range, they obtained relative 
linewidths jq/loq rs 1/3, very consistent with our results 
above. For these modes, the average Xq = 26 (when nor- 
malized as we have done) are again exactly in line with 
the expectations outlined above. 

This extremely strong and abruptly Q dependent EP 
coupling provides new insight into the relationship of 
limits on T c and lattice instability. [35] Earlier study 
had produced arguments both that incipient lattice in- 
stabilites are helpful [33,34] for T c , and that they are un- 
helpful because the low frequency modes are less useful 
for high T c than higher frequencies (of the order of 27rT c ). 
[36] The lack of explicit dependence of A2 9 on /cp indicates 
that increase in T c with doping will arise only from an in- 
crease of to* (the heavy hole mass increases by 25% from 
x = 0toa; = 0.5)or increase of |M 2g | 2 , cither of which 
will further soften the E 2g modes and move the system 
closer to instability. What changes rapidly with doping 
is the coupling strength of the E 2g modes with Q < 2Uf - 
the same total EP coupling strength is concentrated into 
increasingly fewer bond-stretching modes at lower dop- 
ing levels. Thus in this system (where the two dimen- 
sionality of the electron dispersion is central) change of 



the doping level can lead toward (or away from) a lattice 
instability without any direct affect on T c itself. These 
new features involving extremely strong coupling apply 
to Mgi_ x Al x B2, where unexplained structural anomalies 
occur, as well as to Lii_ K BC. 
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FIG. 1. Fermi surfaces ol Lio.7sBC, with the F point at the 
center of the hexagonal zone. There are four cylinders, analo- 
gous to the two cylinders in MgB2 but downfolded due to the 
doubled unit cell. The non-cylindrical surface arises from the 
weakly coupled n bands. 



160 
140 
120 
100 



> 

e so 

3 60 
40 
20 



°r 



















• 














: ' 


• 






















• • < . •-_ 


• 

L • • . 


• : 










t • 






• 










• • 




'•t 












IS 


in 

• • • • ! 
" " * 








: 1 1 ' ' ~ 


• » s ' 1 

1 '■ . • 

- • . 






• 


! i 


• • - 
' . ~, 

- 






• • 











K M 



r a 



L H 




FIG. 2. Calculated phonon dispersion curves for (top) 
semiconducting LiBC and (bottom) strongly hole-doped 
and metallic Lio.7sBC, with i?2 9 -derived modes connected 
by heavy lines. The primary difference is the extremely 
strong renormalization downward (u>q decreases by ~60%) 
for Q < 2kF', the extreme van Hove singularities along 
r - K, M - T, A - L, H - A are apparent. 
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FIG. 3. Phonon DOS for LiBC and Lio.rsBC (top panel) 
where the primary difference is the disappearance of modes 
in the range 135 - 155 meV upon doping (from Fig. 1, they 
mostly lie in the 80-90 meV region). The shape of F(u>) and 
a 2 F(u>) for L10.75BC (bottom panel) revealing the extremely 
strong coupling to phonons in the 85 - 100 meV range. 
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FIG. 4. Characteristics of the model 2D electrons (see text) 
strongly coupled to E2 g modes. Bottom panel: E2 g branch 
versus Q; reference branch Hq and renormalized branch luq as 
in Eqs. (l)-(3), for three different values of 2&f- Top panel: 
Corresponding values of the mode Xq given by Eq. 4, for the 
same three values of 2/cf. 



